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ABSTRACT 

We have continued our long term study of the double-neutron-star binary 
pulsar PSR B 1534-1- 12, using new instrumentation to make very high precision 
measurements at the Arecibo Observatory. We have significantly improved our 
solution for the astrometric, spin, and orbital parameters of the system, as well 
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as for the five "post-Keplerian" orbital parameters that can be used to test grav- 
itation theory. The results are in good agreement with the predictions of general 
relativity. 

With the assumption that general relativity is the correct theory of gravity in 
the classical regime, our measurements allow us to determine the masses of the 
pulsar and its companion neutron star with high accuracy: 1.3332 ± O.OOIOM© 
and 1.3452 ± O.OOIOMq, respectively. The small but significant mass difference 
is difficult to understand in most evolutionary models, as the pulsar is thought 
to have been born first from a more massive progenitor star and then undergone 
a period of mass accretion before the formation of the second neutron star. 

PSR B1534-I-12 has also become a valuable probe of the local interstellar 
medium. We have now measured the pulsar distance to be 1.02 ±0.05 kpc, giving 
a mean electron density along this line of sight of 0.011 cm~^. We continue to 
measure a gradient in the dispersion measure, though the rate of change is now 
slower than in the first years after the pulsar's discovery. 

Subject headings: pulsars: individual (PSR B1534+12) — gravitation — binaries: 
close — stars:distances 



1. Introduction 

The discovery (Hulse & Taylor 1975) of the first pulsar in a binary system, PSR B1913+16, 
introduced the possibility of powerful new experimental tests of gravity in the strong field 
and radiative regimes. Over a quarter century, measurements of this pulsar were shown to 
be in excellent agreement with the predictions of general relativity (Taylor & Weisberg 1989; 
Damour & Taylor 1991; Taylor 1994). A dozen years ago the discovery (Wolszczan 1991) of 
PSR B1534+12 — a bright, relatively nearby pulsar, also with a neutron-star companion — 
promised the independent opportunity to confirm the tests of gravity done earlier with 
PSR B1913+16. Furthermore, B1534+12's favorable, nearly edge-on orbital geometry al- 
lowed new experimental tests that were complementary to those done with PSR B1913-I-16 
(Taylor et al. 1992). 

In earlier work (Stairs ct al. 1998, hereafter Paper 1) we analyzed data taken over a 
six year period with the Arecibo 305 m telescope, the 43-m telescope at Green Bank, West 
Virginia, and the 76-m Lovell Telescope at Jodrell Bank Observatory, U.K. We found that 
those data were in agreement with the predictions of general relativity at the limit of the 
measurement uncertainties: about 1%. We showed that the orbital period was decreasing, 
as expected if the system is losing energy in the form of gravitational radiation. Because 
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the observed orbital decay rate is contaminated by kinematic effects which depend on the 
unknown pulsar distance (Damour & Taylor 1991; Camilo, Thorsett & Kulkarni 1994), the 
radiation predictions of general relativity cannot be tested with this pulsar to better than 
~ 30% without an independent distance measurement. However, with the assumption that 
general relativity is the correct theory of gravity, we were able to invert the kinematic model 
to measure the pulsar distance to 20%. 

Over the past four years, we have continued our studies of PSR B1534+12, using the 
recently upgraded Arecibo telescope together with a new data acquisition system (Stairs 
et al. 2000a) that fully removes the dispersive effects of the interstellar medium from the 
signal. In addition to extending the length of the data set by more than half, these new 
data are of substantially higher quality than those available for Paper I. In this paper, we 
report on timing studies based on the complete 11.5-year Arecibo data set. We present the 
improved tests of general relativity that are now possible, give a pulsar distance accurate to 
5%, and report on dispersion measure variations. In subsequent papers, we will describe the 
analysis of secular variations in the pulsar radiation pattern that have been interpreted as 
evidence for geodetic precession of the pulsar, and we will discuss single-pulse studies. 



2. Observations 

In this paper we discuss results of observations made with two different observing sys- 
tems at the Arecibo telescope. Many of the details have already been discussed in Paper I, 
and are only summarized here and in Table 1. 

Observations were made with the 305-m Arecibo telescope both before its upgrade 
(1990-1994) and after its upgrade (1997-2002). For our purposes, the most important post- 
upgrade change was a significant reduction in the system temperature at 1400 MHz, leading 
to an improvement in the signal-to-noise ratio achieved for a given integration time and 
observing bandwidth. During the period that the Arecibo telescope was unavailable, obser- 
vations were made with the National Radio Astronomy Observatory (NRAO) 43-m telescope 
at Green Bank, and with the 76-m Lovell telescope at Jodrell Bank Observatory, U.K. The 
latter two datasets were used in Paper I, but, with the abundance of high-quality post- 
upgrade Arecibo data, they contribute nothing to the present solution and are therefore 
disregarded in this work. 

Two observing systems were used. The Princeton Mark III system (Stinebring et al. 
1992) is an "incoherent" system, in which the incoming signal is decomposed into spectral 
channels with a filterbank. The effective time resolution is dominated by dispersion smearing 
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of the signal during its propagation through the ionized interstellar medium, with some 
additional instrumental smoothing after detection of the signal. The signals in each frequency 
channel were folded at the instantaneous topocentric pulsar period for three to five minutes, 
then the channels were shifted to account for interstellar dispersion and summed to produce 
a single total-intensity (summed-polarization) pulse profile for each integration. 

For the new post-upgrade data that we introduce here, wc observed with a "coherent" 
system: the Princeton Mark IV instrument (Stairs et al. 2000a). For each sense of circu- 
lar polarization, a 5 MHz bandpass was mixed to baseband using local oscillators in phase 
quadrature. The four resulting signals were low-pass filtered at 2.35 MHz, sampled at 5 MHz 
with 4-bit quantization, and written to a large disk array. Some of the 1400 MHz observations 
were taken with 10 MHz of bandwidth and 2-bit quantization. In the offline processing, the 
undetected signal voltages were dedispersed using the phase-coherent technique described 
by Hankins & Rickett (1975). After amplitude calibrations were applied, self- and cross- 
products of the right- and left-handed complex voltages yielded the Stokes parameters of 
the incoming signal. These products were folded at the topocentric pulsar period using 1024 
phase bins, and pulse TOAs were determined from the resulting total-intensity profiles. A 
summary of the more important parameters and statistics of all combinations of observing 
systems and telescopes is presented in Table 1. Many of the post-upgrade data come from 
observations on one or two days at quasi-regular intervals; in addition we conducted 12-day 
intensive "campaign-style" observations in each of the summers of 1998, 1999, 2000 and 
2001, with full orbital coverage at each epoch. This mixture of data permits a high-quality 
timing solution and allows us to monitor, at each campaign epoch, the pulse shape changes 
due to geodetic precession, which we will discuss in a subsequent paper. 

As in Paper I, we used the same TOA-fitting procedure for all data sets. Each observed 
profile was fitted to a standard template, using a least-squares method in the Fourier trans- 
form domain (Taylor 1992) to measure its time offset. The offset was added to the time of 
the first sample of a period near the middle of the integration, thereby yielding an effective 
pulse arrival time. A different standard template was used for each observing system and 
frequency; they were made by averaging the available profiles over several hours or more. 
Uncertainties in the TOAs were estimated from the least squares procedure, and a minimum 
error, dependent on the observing system and frequency, was required for each TOA. In 
addition, a cut was estabhshed such that all Mark IV TOAs with least-squares uncertain- 
ties greater than 10 /xs (roughly 20% of the data) were discarded; tests show that they do 
not significantly affect the best-fit solution or the parameter uncertainties, and discarding 
these data points allows us to better assess the data quality (§3.3). The observatory's local 
time standard was corrected retroactively to the UTC timescale, using data from the Global 
Positioning System (GPS) satellites. 
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3. Data Analysis 



3.1. The Timing Model 



A pulse received on Earth at topocentric time t is emitted at a time in the comoving 
pulsar frame given by 



Here to is a reference epoch and is the offset between the observatory master clock 
and the reference standard of terrestrial time. The dispersive delay is D/p, where D — 
DM/2.41 X 10-"^, with the dispersion measure DM in cm ^pc, the radio frequency / in 
MHz, and the delay in seconds. The Ar©, A^©, and Asq terms are propagation delays and 
relativistic time adjustments for effects within the solar system, and An, Ae and A5 are 
similar terms accounting for phenomena within the pulsar's orbit. 

In Paper 1, we set out the equations that define the three orbital terms A^, Ae, and 
As in terms of ten parameters. Five are the standard Keplerian parameters: the period Ph, 
projected semi-major axis x = aisini/c, eccentricity e, longitude of periastron uj, and time 
of periastron Tq. The other five are Post Keplerian (PK) parameters: the rate of advance of 
the periastron uj, orbital period derivative Pb, time dilation and gravitational redshift 7, and 
range r and shape s of the Shapiro time delay. These quantities, in conjunction with the 
time-variable dispersion measure, a simple time polynomial to model the spin of the pulsar, 
and astrometric parameters to model the propagation of the signal across the solar system, 
constitute the free parameters to be fit in the theory-independent timing model. 

Within the framework of a particular theory of gravity, the five PK parameters can be 
written as functions of the pulsar and companion star masses, mi and m2, and the well- 
determined Keplerian parameters. In general relativity the equations are as follows (see 
Damour & Deruelle 1986; Taylor & Weisberg 1989; Damour & Taylor 1992): 



T = t-to + Ac- D/p + A^o + Aeq - Asq - Ar 
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Here the masses mi, m2, and M = mi + m2 are expressed in solar units, Tq = GMq/c^ — 
4.925490947 //s, G the Newtonian constant of gravity, and c the speed of hght. 

3.2. Arrival Time Analysis 

We used the standard tempo analysis software (Taylor & Weisberg 1989, see also 
http://pulsar.princeton.edu/tempo) together with the JPL DE200 solar-system ephemeris 
(Standish 1990) to perform a least-squares fit of the measured pulse arrival times to the tim- 
ing model. Results for the astrometric, spin, and dispersion parameters of PSR B1534+12 
are presented in Table 2. The DM was fit in two blocks, pre-and post-upgrade, with different 
base values and time derivatives in each block. Because of systematic errors (see §3.3 below), 
the Mark 111 430 MHz data were only used to fit the pre-upgrade DM and DM derivative, 
and were not incorporated into the rest of the timing solution. This meant that an itera- 
tive approach was required to find the best-fit DMs: first the spin and orbital solution was 
held constant while the DM parameters were fit using the whole dataset, then the Mark III 
430 MHz data were dropped as the best overall timing solution was obtained from the other 
data, holding the pre-upgrade DM parameters fixed. This new timing solution was then used 
to refine both sets of DM parameters, until no significant changes resulted. The DM and 
DM-derivative uncertainties reported in Table 2 are the errors derived while holding the rest 
of the timing solution fixed. Fitting for post-upgrade DM and DM-derivative simultaneously 
with the rest of the timing solution yielded similar values, but somewhat larger uncertainties 
due to covariances with other parameters. Figure 1 shows the time variation of the best-fit 
DM, with the DM and first derivative from Paper I shown for comparison; the pre-upgrade 
values are consistent within the uncertainties. We note that the offset in DM base values 
between the pre- and post-upgrade fits is due mostly to the use of different standard profiles 
with necessarily different frequency alignments; thus this offset should not be interpreted as 
physically meaningful. The change in DM derivative, however, is real, and is discussed in 
§4.2 below. An arbitrary time offset was also allowed between the pre- and post-upgrade 
data sets to accommodate the different standard profiles and changes in the instrumental 
signal path length. 

Two models of the pulsar orbit were used to fit the data. The theory-independent 
"DD" model (Damour & Deruelle 1986) treats all five PK parameters defined in §3.1 as 
free parameters in the fit. The alternate "DDGR" model (Taylor 1987; Taylor & Weisberg 
1989) assumes that general relativity is correct and uses equations 2 through 6 to tie the PK 
parameters to M = mi -|- m2 and m2; consequently it requires only two post-Keplerian free 
parameters. 
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Table 3 presents our best-fit orbital parameters. Uncertainties given in the table are 
approximately twice the formal "1 cr" errors from the fit; we believe them to be conservative 
estimates of the true 68%-confidence uncertainties, including both random and systematic 
effects. For each model, we list the Keplerian parameters and the relevant post-Keplerian 
parameters. We also include, in italic numbers, the computed PK parameters derived from 
the measured masses in the DDGR fit; these computed parameters are in excellent agreement 
with the DD fit values. As in Paper I, the DDGR solution includes the "excess Pf," parameter, 
which accounts for an otherwise unmodeled acceleration resulting from galactic kinematics. 
We note that the measured orbital period derivative, — 0.137 ± 0.003 x 10~^^, is very close 
to the sum of the DDGR model's intrinsic value, —0.1924 x 10~^^, and the fitted excess, 
+0.055 ± 0.003 X 10-^2 

Figure 2 shows the post- fit residuals for the pre- upgrade 1400 MHz data, and the post- 
upgrade 1400 and 430 MHz data, plotted as functions of date. Interstellar scintillation causes 
the TOA uncertainties to vary somewhat, but in the interest of clarity we have not drawn 
error bars on these residual plots. Figure 3 illustrates the averaged post-fit residuals for the 
same datasets, plotted as functions of orbital phase. 

3.3. Comparison of Coherent and Incoherent Observing Systems 

As discussed in Paper I, solutions for the fitted parameters of PSR B1534+12, as de- 
termined from the pre-upgrade 430 MHz data, are biased by small systematic errors in the 
TOAs, likely caused by imperfect post-detection dispersion removal combined with vari- 
able interstellar scintillation patterns. These measurements are unreliable at the several 
microsecond level required for high-precision timing of this pulsar. 

One of our major goals in building the Princeton Mark IV instrument (Stairs et al. 
2000a) was to improve the timing accuracy at low frequencies for fast systems such as 
PSR B1534+12, by minimizing any dispersion-related systematic errors in TOAs. We have 
determined that the new instrument does indeed appear to have circumvented the residual 
dispersion problems that plague the Mark III data at low frequencies. In Figure 4 we follow 
Taylor & Weisberg (1989) and Paper I in examining the statistical properties of the post- 
fit residuals (relative to the solution of Tables 2 and 3) to see if they "integrate down" as 
7i-i/2 when n consecutive values are averaged. The pre-upgrade 1400 MHz data and both 
post-upgrade data sets are fairly consistent with the expected slope, while the pre-upgrade 
430 MHz data deviate significantly. This confirms that we have been justified in leaving the 
pre-upgrade 430 MHz data out of our preferred timing solution, that these Mark III 430 MHz 
data indeed contain systematic errors related to the incoherent dispersion removal, and that 
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the new Mark IV instrument greatly reduces these systematics, as we had hoped. The 
minimum TOA uncertainty requirement mentioned above was introduced to compensate for 
the small remaining deviations from the law visible in the three good datasets. 

4. Discussion 

4.1. Observed Change in Orbital Period: Distance to PSRB1534+12 

The observed P5 is measured in the reference frame of the solar system barycenter, and 
must be corrected to the ccnter-of-mass frame of the binary pulsar system before it can be 
compared to the value predicted by GR. The largest correction is kinematic, involving the 
relative acceleration of the two reference frames. It can be broken down into the vertical 
acceleration in the Galactic potential, acceleration in the plane of the Galaxy, and an ap- 
parent centripetal acceleration due to the transverse velocity of the pulsar binary (Damour 
& Taylor 1991): 



Here ttz is the vertical component of galactic acceleration, / and b the pulsar's galactic 
coordinates, i?o and vq the Sun's galactocentric distance and galactic circular velocity, /x and 
d the pulsar's proper motion and distance, and P = d/Ro — cosl. As we do not have a precise 
distance to the pulsar through either a timing or an interferometric parallax measurement, 
the best independent distance estimate still comes from the Taylor & Cordes (1993) model 
of the free electron content of the Galaxy; this model puts the pulsar at a distance of 
0.7±0.2kpc. At this distance the Kuijken & Gilmore (1989) model of the Galactic potential 
yields an estimate of Qz/c— (1.60 ±0.13) x lO^^^s"^. We assume Vo = 222 ±20 km s~^ and 
Ro = 7.7±0.7kpc, as in Damour & Taylor (1991). Then, summing the terms in equation (7) 
and multiplying by Pb, we find the total kinematic bias to be 



The uncertainty in this correction is dominated by the uncertainty in distance, which is only 
roughly estimated by the Taylor and Cordes model. The shght decrease in the uncertainty 
from that given in Paper I results from an increase by a factor of ten in the precision of the 
proper motion measurement. 

Our measurement of the intrinsic rate of orbital period decay is therefore 
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The uncertainty is completely dominated by the uncertainty on the kinematic correction, 
which is nearly a factor of 4 larger than the measurement uncertainty on P^°^^. In GR, 
the orbital period decay due to gravitational radiation damping, (Pfe)*^^, can be predicted 
from the masses mi and m2 (eq. 4), which in turn can be deduced from the high precision 
measurements of cu and 7. As listed in Table 3, the expected value is 

/ . N GR 

(pA = -0.192 X 10-^^ (10) 

As in Paper I, our measured value differs from this prediction, now by some 1.7 standard 
deviations, and again, assuming that GR is the correct theory of gravity, we can derive 
the true distance to the pulsar from the "excess Pb" parameter and equation 7 above (Bell 
& Bailes 1996). Our improved distance is li = 1.02 ± 0.05 kpc (68% confidence limit). 
The uncertainty is still dominated by the measurement uncertainty of {Pb)°^^, rather than 
uncertainties in the galactic rotation parameters or the acceleration a^, though the Galactic 
model uncertainties will ultimately limit the distance measurement to a few percent accuracy. 
The timing parallax for this system is still not significantly measured, but is constrained to 
be less than 1.5 mas, in good agreement with the GR-derived result. 

Our new distance result is consistent with the 1.1 ± 0.2 kpc determined in Paper I, 
but we have improved the uncertainty by a factor of 4. Our previous distance led to a 
downward revision of the estimated double-neutron-star inspiral rate visible to gravitational- 
wave observatories such as LIGO. The new measurement leads to a small (15%) increase 
in the number density of similar systems in the local universe, with the uncertainty now 
dominated by the uncertain scale height of such binaries (e.g., Kalogera et al. 2001). 

4.2. Dispersion Measurements and the Local Interstellar Medium 

As noted in §3.2, we find a significantly different DM and DM derivative in the post- 
upgrade era than before the upgrade (Figure 1). We argue that while the change in reference 
DM value is largely due to the use of new standard profiles with a slightly different frequency 
alignment, the difference in the rate of change (more than a factor of 4 slower since the 
upgrade) is physical. It is true that the profile at both frequencies is undergoing a secular 
evolution due to geodetic precession of the pulsar's spin axis (Arzoumanian 1995, Stairs 
et al., in prep.) and thus one might suspect a contribution to the DM derivative from 
different temporal evolution at the two observing frequencies. We explore this possibility by 
calculating the DM difference between our 1999 May and 2001 June observing campaigns, 
for the best-fit post-upgrade DM derivative versus the best-fit pre-upgrade DM derivative, 
finding that if the pre-upgrade DM trend had continued, the 2001 June DM would have been 
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smaller by about 0.0005 cm~^ pc. This large a difference would imply a difference in offset 
between the 430 MHz and 1400 MHz profiles of 10 fis, or 0.27 bins in a 1024-bin profile. By 
comparing the actual cumulative profiles at these epochs using the same cross-correlation 
algorithm used to calculate TOAs, we estimate that the 430-1400 MHz offset has shifted by 
no more than about 0.05 bins. As the TOAs arc determined by the strong main peak of the 
pulse and the most noticeable evolution is taking place in the low-level emission near the base 
of the pulse (Stairs et al. 2000b) this is perhaps not surprising. We therefore conclude that 
the observed changes in DM slope do represent changes in the structure of the interstellar 
medium between the Earth and the pulsar. 

With our current data set, it is not possible to construct a satisfactory structure func- 
tion with which to characterize the length scales of turbulence in the interstellar medium 
(e.g., Phillips & Wolszczan 1991). Because of strong refractive scintillation at 1400 MHz, the 
high-frequency data set remains quite sparse (Figure 2). Even though we may use a wider- 
bandwidth instrument such as the recently commissioned 100-MHz "WAPP" spectrometer 
for future 1400 MHz observations, the broadband nature of the refractive interstellar scin- 
tillation may make it difficult to accurately sample the dispersion measure variations with 
resolution finer than roughly 6 months. 

4.3. The Neutron Stair Masses: Testing Evolutionary Theory 

The best mass estimates for the two neutron stars in the PSR B1534-I-12 system come 
from the DDGR timing model. The measured values are mi = 1.3332 ± 0.0010 M0 and 
7722 = 1.3452 ± 0.0010 Mq. The masses derived from the DD model and equations 5 and 6 
are in good agreement with these values, though the measurement uncertainties are much 
larger (Figure 5). 

Although the masses of the two neutron stars are very similar, it is now clear that the 
pulsar is significantly less massive than its companion. This is contrary to initial expec- 
tations from binary evolution. The pulsar is rapidly spinning, and has the low magnetic 
field characteristic of a "recycled pulsar" that has been spun-up by mass transfer from a 
companion star. Assuming a monotonic relation between progenitor mass and neutron star 
mass, we therefore might predict that the pulsar would be more massive at birth, and that 
the mass difference would only increase during the later evolution of the system. 

This perhaps naive prediction has already been challenged by the binary pulsar system 
PSR B2303-I-46. Recently, the companion in this close, eccentric binary was shown to be 
a white dwarf rather than a second neutron star (van Kerkwijk & Kulkarni 1999). The 
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eccentric orbit implies that the white dwarf formed first, from the star that was originally 
more massive. Presumably, mass transfer to the initially less massive star pushed it above 
the minimum mass needed to form a supernova and a neutron star. 

A similar process may have been important in the formation of the PSR B 1534+ 12 
system. Although the pulsar was formed from the initially more massive star, mass transfer 
to its companion star probably resulted in a mass inversion before either neutron star was 
formed. 

To move beyond the simplest (or most naive) predictions about the relative masses of 
the pulsar and companion in the B1534+12 system requires significant improvements in our 
theoretical understanding of several areas of stellar and binary evolution. First, we must 
understand the relationship between progenitor and neutron star masses, which may not be 
purely monotonic. Then we must understand how the mass lost by the pulsar's progenitor 
or gained by the companion's progenitor affected the evolution of the stellar cores and the 
amount of fallback during the supernova events. Finally, we must understand in some detail 
the mass transfer that spun up the pulsar. Although these are all challenging problems, there 
is some hope for future guidance from the growing number of precisely measured neutron 
star masses. 



4.4. Test of Relativity 

PSR B1534+12 permits the second test of general relativity based on the to, 7, and Pb 
parameters of a binary pulsar system, and adds significant measurements of the Shapiro-delay 
parameters r and s. The left-hand sides of equations (2-6) represent quantities measured in 
a theory-independent fashion and listed in the DD column of Table 3. If GR is consistent 
with these measurements and there are no significant unmodeled effects, the five curves 
corresponding to equations (2-6) should intersect at a single point in the mi-m2 plane. 
These curves arc presented in Figure 6, in which a pair of lines delimit the 68% confidence 
limit for each PK parameter (a single line is drawn for uj, whose uncertainty is too small to 
show). It is clear that the ci;, 7, r, and s curves intersect, though r is still poorly measured. 
The curve obtained from the observed DD value of can be made to intersect the others, as 
discussed above, by setting the pulsar distance to 1.02±0.05 kpc rather than the 0.7±0.2 kpc 
estimated from the dispersion measure. A filled circle at mi = 1.3332 MQ,m2 — 1.3452 Mq 
marks the DDGR solution of Table 3, and its location on the 00 line agrees to within 0.05% 
with the measured DD values of 7 and s. This provides a highly precise test of the validity 
of general relativity using only non-radiative timing parameters, an important complement 
to the mixed uj^-Pi, test provided by PSR B1913-I-16. 
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Fig. 1. — The dispersion measure (DM) used in the fits presented in Tables 2 and 3. Different 
values and derivatives are needed for the pre- and post-upgrade data sets. The DM and its 
variation obtained in Paper I are shown for comparison; they match our current results 
within the uncertainties. The points with error bars illustrate the quality of DM fits over 
short intervals: the horizontal errors bars indicate the domain of each fit, while the vertical 
error bars represent the uncertainties derived in each time bin. 
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Fig. 2. — Post-fit residuals versus date for (a) the pre-upgrade (Mark III) 1400 MHz data, (b) 
the post-upgrade (Mark IV) 1400 MHz data, and (c) the post-upgrade (Mark IV) 430 MHz 
data. 
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Fig. 3. — Average post-fit residuals as a function of orbital phase for (a) the pre-upgrade 
(Mark III) 1400 MHz data, (b) the post-upgrade (Mark IV) 1400 MHz data, and (c) the 
post-upgrade (Mark IV) 430 MHz data. 
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Fig. 4. — Root-mean-square residual versus number of consecutive residuals averaged, for 
the four data sets: prc-upgradc (Mark 111) and post-upgrade (Mark IV) data at each of 1400 
and 430 MHz. The short-dashed line indicates the expected slope of —1/2 for uncorrelated 
residuals. All data sets except the pre-upgrade 430 MHz data are in reasonable agreement 
with this prediction. 
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Fig. 5. — The Shapiro delay parameters. Contours indicate the 68% and 95% confidence 
levels for the companion mass and the cosine of the orbital inclination. The round dot 
represents the best-fit DD solution in Table 3 while the triangle represents the optimal 
DDGR solution, which is well within the 68% contour of the DD solution. 
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Fig. 6. — Mass-mass diagram for the PSR B1534-I-12 system. Labeled curves illustrate 
68% confidence ranges of the DD parameters hsted in Table 3. The filled circle denotes the 
component masses according to the DDGR solution. A kinematic correction for assumed 
distance d = 0.7 ± 0.2 kpc has been subtracted from the observed value of Pb, the uncer- 
tainty on this kinematic correction dominates the uncertainty of this curve. A slightly larger 
distance removes the small apparent discrepancy between the observed and predicted values 
of A. 



Table 1. Parameters of the four observing systems. 



Mark III Mark III Mark IV Mark IV 



Frequency (MHz) 1400 430 1400 430 

Bandwidth (MHz) 40 8 5 5 

Spectral Channels 32 32 1 1 

Dedispersing system . . . incoherent incoherent coherent coherent 

Time resolution (//s) ... 97 329 38 38 

Integration time (s) ... 300 180 190 190 

Dates 1990.8-94.1 1990.7-94.2 1998.6-2002.3 1997.9-2002.3 

Number of TOAs 1185 2311 639 1635 

Median cttoa (/^s) 6.1 6.0 6.7 4.2 
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Table 2. Astrometric, spin, and dispersion parameters for PSR B1534+12? 



Parameter 


Value 


Rigiit ascension, a (J2000) 


15^37°'09!960312(10) 


Declination, S (J2000) 


11° 55' 55'.'5543(2) 


Proper motion in R.A., /la (masyr~^) . 


1.32(3) 


Proper motion in Dec, /is (masyr^^) . . 


-25.12(5) 


Parallax, tt (mas) 


< 1.5 


Pulse period, P (ms) 


37.9044407982695(4) 


Period derivative, P (10~^®) 


2.422622(3) 


Epoch (MJD) 


50261.0 


Dispersion measure, DM, 1 (cm~^pc) . . 


11.61944(2) 


DM derivative 1 (cm~^pcyr~^) 


-0.000316(10) 


DM Epoch 1 (MJD) 


48778.0 


DM 2 (cm-3pc) 


11.61634(3) 


DM derivative 2 (cm~^pcyr~^) 


-0.000073(7) 


DM Epoch 2 (MJD) 


51585.0 


Galactic longitude / (deg) 


19.8 


Galactic latitude b (deg) 


48.3 


Composite proper motion, /i (masyr~^) 


25.15(5) 


Galactic position angle of ji (deg) 


239.16(8) 



^Figures in parentheses are uncertainties in the last digits quoted, 
and italic numbers represent derived quantities. 
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Table 3. Orbital parameters of PSR B1534+12 in the DD and DDGR modelsf 





DD model 


DDGR model 


Orbital period, Ph (d) 


0.420737299122(10) 


0.420737299123(10) 


Projected semi-major axis, a; (s) 


3.729464(2) 


3.7294641(4) 


Eccentricity, e 


0.2736775(3) 


0.27367740(14) 


Longitude of periastron, ui (deg) 


274.57679(5) 


274.57680(4) 


Epoch of periastron, Tq (MJD) 


50260.92493075(4) 


50260.92493075(4) 


Advance of periastron, u (degyr"-*-) 


1.755789(9) 


1.7557896 


Gravitational redshift, 7 (ms) 


2.070(2) 


2.069 


Orbital period derivative, {Pb)°^^ (10~^^) 


-0.137(3) 


-0.1924 


Shape of Shapiro delay, s 


0.975(7) 


0.9751 


Range of Shapiro delay, r (/xs) 


6.7(1.0) 


6.626 


Derivative of x, \x\ (10~^^) 


< 0.68 


< 0.015 


Derivative of e, e (10~^^s~^) 


< 3 


< 3 


Total mass, M — rrii + m2 (Mq) 




2.678428(18) 


Companion mass, m2 (Mq) 




1.3452(10) 


Excess Pb (10-12-) 




0.055(3) 



^Figures in parentheses are uncertainties in the last digits quoted. Italic numbers 
represent derived parameters, assuming general relativity. 



